A model of the T-dependent pseudogap and its competition with superconductivity in 

copper oxides. 
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Results for pseudo-gaps are obtained from a band model, where the stability of the gap depends on 
the amplitudes of vibrational displacements, or magnetic moments, and their coupling to electrons. 
A one-particle gap is favored by normal thermal excitations of phonons or spin waves. Another 
gap can be generated by spontaneous waves at lower temperature, if the electronic energy gain 
overcomes the elastic/magnetic energy needed for increased amplitudes of the oscillations. This 
state is characterized by charge or spin density waves. The pseudogap has many features in common 
with the superconducting gap, and the model lends support to the interpretation that the pseudogap 
is a precursor of, and competes with, superconducting pairing. 

PACS numbers: 74.20.-z,74.20.Fg,74.20.Pq 



The understanding of high-T^ cuprates has increased 
a lot since the discovery of these systems more than 20 
years ago although a convincing mechanism behind 
the high superconducting Tc is not found. A pseudogap 
with transition temperature T* > Tq at low hole doping 
seems to be characteristic for all cuprates [2-4] . Ab-initio 
band calculations for long supercells of copper oxides con- 
taining phonon distortions or spin waves show a partial 
gap (pseudogap) in the density-of-state (DOS) 0,|g|. The 
energy position of this gap depends on the length of the 
supercell, i.e. on the wave length of the periodic po- 
tential perturbation induced by the phonon/spin wave. 
Hence, if the one-particle gaps appear at the Fermi en- 
ergy, Ep, then there must be a correlation between dop- 
ing (determines Ep) and wave length. The appearance of 
the gap can be understood from the nearly-free-electron 
model (NFE), which has also been used to interpolate 
ab-initio results Q • The pseudogap is an important part 
of cuprate physics, but the question if it helps or com- 
petes with superconductivity is not settled. On the other 
hand, it can be argued that artificial periodic waves such 
as given by periodic doping or oxygen orderingwill lead 
to higher Tc , if the waves are tuned correctly [8|, Q . 

Presented here is a many-particle extension of the 
NFE band model, which can be used for the simulation 
of the T-dependence of the pseudogap. The model is 
based on one-dimensional potential perturbations from 
phonons or anti-ferromagnetic (AFM) spin waves. The 
criterion for having gaps at Ep is made from the esti- 
mated total energy (kinetic electronic energy and energy 
of phonon/spin fluctuations) for an interacting electron- 
phonon (electron-spin) system, with many-body coupling 
parameters A or A s / for phonons or spinfluctuations, re- 
spectively. A continuous metallic band can be gapped 
because of a periodic potential perturbation, as in the 
appearance of a gap for semiconductors [1(| . The poten- 
tial perturbation is such that two regions of k-space, at 
k and -k, are affected equally. Only waves which mod- 
ify electronic states around the Fermi surface (FS) are of 
interest, since changes of the bands far from the Fermi 
energy (Ep) make no change in the total kinetic energy. 



This puts a constraint on the g-vectors of the perturba- 
tions. The results of the model imply similar mechanisms 
for the pseudogap and the superconducting gap in high- 
Tc cuprates. 

Phonons and spin fluctuations are normally ex- 
cited thermally following the Bose-Einstein occupation, 
g(T,u), of the phonon- or spin wave density-of-states 
(DOS), F(oj) or F m (uj), respectively. The averaged 
atomic displacement amplitude for phonons, u, can be 
calculated as function of T [l(| [H|. Approximate re- 
sults make u 2 z — > 'ihwo/^K at low T ("zero-point mo- 
tion", ZPM) and u\ -> 3k B T/K at high T ("thermal 
excitations"), where ujp> is a weighted average of F(ui). 
The force constant, K = Maui 2 , where Ma is an atomic 
mass, can be calculated as K = d 2 E/du 2 (E is the to- 
tal energy). The corresponding relations for averaged 
fluctuation amplitudes of the magnetic moments, m, are 
the same, but without the polarization factor 3 and with 
the replacement of K with K m = d 2 E/dm 2 [12]. The 
force constants K and K m do not change with the os- 
cillation amplitudes for harmonic oscillations. The time 
scales of phonons and electrons are sufficiently different 
for adiabatic relaxation of electrons and band gaps. The 
energy of an atomic oscillation, U, has an elastic contri- 
bution because of u, and a kinetic contribution because 
of the velocity, v. The time dependence of the sum, 
217 (t) = Ku 2 cos 2 (ujt) + M av 2 sin 2 (wt) , is a constant in 
the harmonic approximation. J7(0) = ^Ku 2 , where u 
refers to the maximal atomic displacement, permits us to 
calculate U from static conditions, and to identify Ku 2 as 
an ingredient of the standard expression for A. Acoustic 
phonons are more efficient for having a clear gap, while 
optical phonons sometimes have u=0 everywhere and a 
smeared time average of the band gap. 

Thermally activated phonons and spin waves create 
some disorder and will change the electronic states and 
the DOS, N(E) [l^]. Individual phonon- or spin- waves 
which cause a gap close to Ep in the normal state are par- 
ticurarly interesting for high-Tp copper oxides p|. The 
gap is largest in the CuO bond directions, while the FS 
remains sharp in the diagonal direction, to form an FS- 



2 



arc |l4[. Band calculations show that atomic distortions 
of phonons create a periodic potential along a chain of 
atoms, which can be modeled by 



V(x) 



V q e- lx - q 



(1) 



for phonon propagation along x with wave vector q (5l.[To|. 
The maximum amplitude, V q , can be obtained from ab- 
initio band results as half of the band gap, or as the max- 
imum difference of the potential shift within the unit cell. 
A spin wave makes an analogous perturbation within the 
spin polarized part of the potential, where the densities 
of opposite spins differ by a phase factor of tt. The result 
is an AFM spin configuration with wave length given by 
2ir/q. Phonons and spin waves can be considered sepa- 
rately, but an important spin-phonon coupling (SPC) in 
the cuprates leads to unusual properties [g, Il4l |15| . The 
following description, based on phonon excitations, can 
be adopted to spin waves and SPC. 

The NFE band has a gap E g — 2V q at some "zone- 
boundary" k = G/2 in one dimension (1-D) (EOS- (All 
q-vectors are concerned initially, but a few of them are 
particular because of energy gains from the gap at Ep.) 
The general band dispersion as function of k 



e = -(k + (k - G) ± ^(k 2 - (k - G) 2 ) 2 + AV 2 ) (2) 

can be made simpler near Ep where e ~ const. ■ k, when 
the k-point, k, is measured from the zone boundary. 



£ = ± 



(3) 



if e is much smaller than the band width, W. The wave 
function at k — is commensurate with V(x), while for 
other k-points, when e is below Ep, it is not. 

The normal free electron dispersion, e = e, is recovered 
for V q — 0, and N(e) = N/\de/de\, becomes constant and 
equal to N. The gapped N is zero for e = Ep ±V q . At 
very low T there is practically no thermal occupation of 
phonons. The Fermi-Dirac occupation / is essentially a 
step function at Ep. 

The propensity for having small or large gaps can 
be found through ab-initio band calculations, since one- 
particle bands are more or less sensitive to the potential 
perturbation. It depends on the material if this is go- 
ing to be an important effect or not. The total energy 
is reduced if a gap is formed at Ep because of the per- 
turbation, but it is not lower than the total energy for 
the system without phonon distortion. In the following 
it is argued that an additional effect, beyond the one- 
particle band mechanism, can lead to spontaneous exci- 
tations of charge/spin wave gaps in systems where A or 
A s / are large. The A's are results of many-body inter- 
actions across Ep, and they are not solely determined 
by the one-particle band structure. The total energy for 
the vibrating system can then be lower than for a static 
system. 



The electron-phonon coupling A is active for energies 
±hw around Ep, where it can be written NM 2 /K [10l | . 
The matrix element M for energies inside the interval 
±hu> can be evaluated as (** (Ep , r)dV (r) / du^J (Ep , r)) , 
which is the first order change in energy caused by the 
perturbation dV(r) for du — > 0. For a finite value of 
u the change in energy will be finite and equal to the 
gap V q , since V q /u is constant for harmonic vibrations. 
Thus, instead of calculating M as a matrix element it is 
convenient to take the value directly from the band gap, 
and M can be written V q /u for energies close to Ep, 
which makes A = NV q 2 /Ku 2 . 

Totally there is a gain in energy if \U\ < \E\. The sys- 
tem will spontaneously increase u of the vibrations (from 
the normal value given by thermal excitations) in such a 
case. Other effects such as electron-electron correlation 
and potential terms add to the energy costs and can pre- 
vent a gap in many systems. With the DOS, U and E 
per unit cell the condition \U\ = \E\ is written 



-Ku 2 = 
2 



e(N(e) - N(e))de 



(4) 



—huj 



for T — 0, where N(e) is the DOS with the gap and 
N(e), the DOS of the normal state, is assumed constant 
within fku around Ep. The integration is to hw since A 
is zero for energies larger than ±Huj. 

With a substitution e 2 = e 2 + V 2 we obtain TV = 



N\e\/Je 2 - V 2 and, 



-Ku 2 = 



Nede- 
huj J —huj 



Ne 2 /^J(e 2 - V 2 )de (5) 



The result is 



and 



Ku 2 = NV 2 ln(2huj/Vq) 



V q = 2hue- 1/x 



(6) 



(7) 



since NV 2 / Ku 2 turns out to be equal to A. The cou- 
pling determines the gap through constant ratios of V q /u. 
This result is similar to the BCS equation for the super- 
conducting gap 0, [I?]]. The present state is different 
from a normal gapped state obtained via standard band 
calculations. It only appears because of electronic inter- 
action near Ep, and it disappears for small A. 

Eqn. [5] has no solutions for Huj < V q , so the concurrent 
state with phonon softening, where all electronic energy 
goes into a renormalization of the phonon with K, uj — > 0, 
is not described. Neither using W instead of hui in eqn [5] 
gives a proper description of a static case, because states 
far below Ep do not have the correct wavelength and 
cannot contribute to V q . 

The model for finding the maximal temperature for 
having a gap (T*) is obtained from eq. UJ but with the 
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FIG. 1: The full line shows 1/A as a function of feT as ob- 
tained from the BCS result for fuo=Q.l eV. The thin and 
broken lines show the corresponding result from eqn. [S] for 
5 different V q 's (0.5, 1, 2, 4 and 8 meV). The result for the 
smallest V q is indistinguishable from the BCS result over this 
temperature range. 

Fermi-Dirac function, /, as the T-dependent weight fac- 
tor for N(e) and N(e), and with the integration in the 
interval [— fraj,huj}. At T* it is required that V q — > 0. 
This is solved numerically through integrations I; 

Ku 2 /NV q 2 = 1/A = I(V q ) = 




l^{e 2 ~V 2 )fde)/V 2 (8) 



for V q — > 0. The ' in the second integral means that 
the energies where |e| < V q are excluded. 

The results, shown in Fig. 1, has the same form as the 
solution for superconductivity: 

k B T* = l.l3fiue- 1/x (9) 

For example, it can be verified from this formula and 
Fig [1] that a A of 0.5 makes fcgT* ss 15 meV when fruj is 
100 meV. 

A non-constant DOS, with N(e) from eqn. Q] inserted 
in eqn. and tkJ being a large fraction W, tends to 
decrease T* . 

For AFM spin waves there is a cost in magnetic en- 
ergy, which in the harmonic approximation can be writ- 
ten U m = ^K m m 2 . The change in potential on some 
site, V m , is positive for one spin and negative for the 
other spin direction. This defines a X s t — NV^/ Km 2 as 
a coupling constant for spin- fluctuations (l2j |. The rest 
of the equations are applicable with A s / replacing A and 
with hui s f being the energy of the spin wave. 

Typical atomic displacements and magnetic moments 
from phonons and spin waves in the normal state can 
be determined from the effective force constants K and 



K m and the gap values, u = \f{N '/ KX)V q and m = 
y/(N/K m \ s f)V q , respectively. 

The charge/spin wave state has low resistivity if the 
gap is complete over the FS, i.e. N(Ef) = 0. An elec- 
tric field E x , applied during a time r, accelerates free 
electrons until their velocity x becomes eE x T/m, where 
e and m are the electron charge and mass, respectively. 
The velocity v x (k) changes from Hk/m to hk/m+x (k 
is the momentum along x), and the free electron band 
from h 2 k 2 /2m to h 2 k 2 /2m+ xk, as shown in fig. [2] 
For a metal, with Ep indicated by the horizontal line 
in fig. [21 occupied states near k = +1 (fci?) are pushed 
upward, while the states near k = —1 (— kp) decrease 
their energy. The net velocity J v x (k)dk ss 2vpx and the 
conductivity is approximately e 2 rNvp /m, as from the 
Boltzmann equation [lfj, |l8j. But the total kinetic en- 
ergy E t = J N(e)ede is enhanced by the field, so the net 
current will decay through resistive dissipation when the 
field is switched off. A band insulator, with a gap at the 
zone boundary, will also have a higher E t from an applied 
field, and by switching off the field, E t decays back to its 
low value for the ground state. 

However, if the free electron band without field (thin 
line in fig. is shifted horizontally by the electric field to 
a new position (bold line) , then E t will remain at the min- 
imum, no decay is possible, so the state with the net cur- 
rent can persist. It is equivalent to moving the origin of k 
to —2m/hx. The electron occupation of the gapped band 
will not be asymmetric as the bold broken line in fig. [21 
but will be symmetric relative to the new origin as for the 
unshifted band. This can be achieved through phonons 
with slightly different momentum so that phonons and 
electrons have a common drift. (Two phonons with 
momentum ±P have energy hflp = K / M A\sin(P a)\ 
and velocity ±ay/ K/Ma\cos(Pcl) \ , where a is the real 
space periodicity. Together they have energy 2Mlp and 
zero velocity. If the phonons change to P + p and 
— P + p (and p << P), there will be a net velocity, 
V p « 2oa/ K / M A\sin(P a)\pa —x. The phonon energy, 
2hilpcos(pa), will not increase, since cos(pa) < 1.) 

A common drift of electrons and phonons implies dif- 
ferent lattice conduction along opposite directions rela- 
tive to an imposed current, which also will be different for 
a state with no current. Effects on lattice and spin- wave 
dependent properties should be observable although very 
small. 

From the selective q-dependence it can be shown that a 
weak magnetic field will destroy the pseudogap. As men- 
tioned, one q generates the gap on the paramagnetic FS 
at Ep. A magnetic field, H, will split the FS into two, one 
for each spin ("up" or "down"), and E F ' d = Ep ± hbH- 
Two independent phonons are required for having the 
gaps optimally on two FS. One has a potential perturba- 
tion at e~ l w~ s >' s , and the other at e - l (?+ s )- x ; where S is 
determined by the band dispersion and (J-bH. The sum 
of these two potentials is 2cos{5 ■ x)e~ lq ' x , and therefore, 
even if there is a modulation given by the cosine function, 
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FIG. 2: Thin, full and broken lines; free electron band without 
and with a gap at the zone limits (vertical lines), respectively. 
Bold full and broken lines: The bands modified by an electric 
field. 

the effective q remains the same and cannot fit optimal 
values for two FS. The resulting gaps do not appear at 
Ep on the two spin-split bands, which will reduce the 
gain in energy. The energy difference, D(H,T), between 
the kinetic energy for the gapped DOS with and without 
field, which is calculated as 

eN(e)(f(e+H, T)+f(e-H, T)-2/(e, T))d, 

-Huj' 

(10) 



increases quadratically with the field amplitude H < V q 
for low and high T (« V q ). This is because the thermal 
occupation can be made more efficiently if Ep is closer to 
the DOS peak (on N) above the gap for "majority" and 
closer to the DOS peak below the gap in the "minority" 
states, than if Ep is in the middle of the gap. 

In conclusion, three mechanisms promote gaps at Ep. 
First, thermal normal state excitations of phonons/spin 
waves generate potential perturbations and weak band 
gaps near Ep. The total energy for the perturbed state 
is not lower than for the static state, but phonons can 
be softened. Second, systems with large A can sponta- 
neously generate phonons/spin waves because of a gain 
in total energy, as has been described above for NFE 
bands. The third mechanism is through superconduct- 
ing pairing. The equations for T* and Tq are similar. 
Also conductivity and the effect of a magnetic field be- 
have in similar ways [T3, EH . Calculations show increased 
A for low doping, which would agree with the evolution 
of T* . But the superconducting Tc, which is maximal 
at a doping of 0.12-0.15 holes/Cu, does not follow the 
doping dependence of A. The shape of the FS and dif- 
ferences in the k, fc'-dependence of A (nesting etc.) are 
probably important for which state will win. In 3-D it 
is hard to imagine that a multitude of density waves can 
coexist to make a full gap everywhere on the FS sphere. 
For 2-D cuprate bands it is simpler since A is only large 
for k-points near X and Y, where a pseudogap can ap- 
pear. These results support the picture that the pseudo- 
gap competes with superconductivity, and that it seems 
to depend on preformed pairs between k and — k states. 
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